Dynamical Transition in Sliding Charge-density Waves with Quenched Disorder 
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We have studied numerically the dynamics of sliding charge-density waves (CDWs) in the presence 
of impurities in d=l,2. The model considered exhibits a first order dynamical transition at a critical 
driving force Fc between "rough" (disorder dominated) {F < F^) and "flat" {F > Fc) sliding phases 
where disorder is washed out by the external drive. The effective model for the sliding CDWs in 
the presence of impurities can be mapped onto that of a magnetic flux line pinned by columnar 
defects and tilted by an applied field. The dynamical transition of sliding CDWs corresponds to the 
transverse Meissner effect of the tilted flux line. 
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I. INTRODUCTION 

The study of the dynamics of an ordered medium 
driven by an external force through quenched impurities 
is relevant to many physical systems. Examples include 
charge-density waves (CDWs) in anisotropic metals j|] 
and flux-line lattices (FLL) in the mixed state of type- 
II superconductors ||]. It has been argued that in both 
these systems impurities destroy the equilibrium long- 
range translational order (LRTO) and pin the medium. 
A driving force F, originating from an electric field or a 
current, can overcome the constraining forces from the 
impurities and cause the medium to slide. At T = 0, 
classical CDWs exhibit a depinning transition at a criti- 
cal driving force Ft from a pinned {F < Ft) to a sliding 
{F > Ft) state. This transition has been described as a 
dynamic critical phenomenon. The nonlinear dynamics 
of the system near the critical point has been studied ex- 
tensively both by numerical simulation and by 4 — e 
expansions Q and is fairly well understood. 

In contrast, the dynamics of driven disordered media at 
large driving force, well above the depinning transition, 
has only recently begun to receive some attention. In the 
sliding state the pinning by impurities is less effective and 
it has been suggested that the medium may recover the 
LRTO at sufficiently large velocity. Recent experiments 
in YBCO as well as simulations of two-dimensional flux 
lattices have indeed shown that the flux array orders at 
large drives. Koshelev and Vinokur Q have described 
this phenomenon as a true phase transition from a flow- 
ing liquid to a moving solid. Whether the ordering of 
driven flux lattices is a true dynamical phase transition 
or a crossover is still an open question. 

Some of us recently addressed this class of questions 
by focusing on a model of driven CDWs in the presence 
of disorder and thermal noise 1^ . The model considered 
in Ref. |^ allows for disorder-induced phase slips of the 
CDW. In three dimensions it yields a dynamical phase 
transition of the sliding CDW from a disordered phase 
with plastic flow to a temporally periodic "moving solid" 



phase with quasi-long-range translational order. In two 
dimensions the moving solid phase is unstable due to the 
proliferation of phase slips. 

In the present paper we focus on a related model 
of the dynamics of sliding CDWs that incorporates 
a nonequilibrium nonlinear term of the Kardar-Parisi- 
Zhang (KPZ) form §], that was neglected in Ref. §]. To 
render the analysis tractable with this addition, we have, 
however, neglected phase slips, and, for most of what fol- 
lows, thermal fluctuations. While the neglect of thermal 
noise has very little effect upon our results, the omission 
of phase slips appears more severe. Indeed, the model 
considered here is by definition an elastic continuum at 
all driving forces and a "liquid phase" with dislocations 
cannot occur. We expect that at sufficiently low tem- 
perature and large driving force the omitted defects will 
make only minor changes to our results, at least on ex- 
perimentally observable time scales. The true (infinite 
time) asymptotics is, however, likely to be affected by 
this omission. 

Interestingly, despite its strong topological constraints, 
our model nevertheless exhibits a phase transition at a 
critical driving force Fc from an isotropic rough (disor- 
der dominated) flowing phase at small driving forces, to 
an anisotropic smooth flowing phase where nonlinearities 
are washed out by the external drive. The average CDW 
velocity changes sharply at the transition, which is ar- 
gued to be first order. The two phases are characterized 
by different values of the roughening exponent governing 
the growth of spatial fluctuations of the CDW phase with 
the size of the system in the direction of the exter- 
nal drive (a;— direction). In the rough phase the phase 
fluctuations grow linearly with L^, 'w{Lx, L±) ~ L^, in- 
dicating that the elastic model breaks down. This sug- 
gests that in a corresponding model that allows for phase 
slips, this phase would be a flowing liquid. Above Fc 
the disorder is washed out in the direction of motion. 
Fluctuations are strongly suppressed in this direction and 
w{L^,L±) ~ L^^ for » L±, with z « 0.85 ± 0.05. 

Our results have two implications for CDW experi- 
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ments at high velocities, provided phase sHp effects are 
sufficiently suppressed at low temperatures. First, above 
the critical force {F > Fc), translational correlations are 
expected to be highly anisotropic, decaying much more 
rapidly transverse to the motion than along it. This im- 
plies a substantial increase in the ratio of widths of Bragg 
scattering peaks. 



d fAk± 
dF \~Kh> 



>0, 



(1.1) 



where Afcj^ and Afcy are peak widths perpendicular and 
parallel to the CDW wave vector, respectively. Secondly, 
within our model, the non-equilibrium ordering transi- 
tion is characterized by a jump discontinuity in the dif- 
ferential conductance Gas = dl/dV. As discussed in 
further detail in the next section, we are able to predict 
only the singular behavior of Gdig, and not the full non- 
linear form of the I(V) curve. 

We expect these two results to survive in varying de- 
grees in models including phase slips (and hence exper- 
iments). The strong anisotropy of the translational cor- 
relations in the high- velocity phase should remain at low 
temperatures. The resistance singularity is expected to 
be more sensitive to phase slips (as, indeed, is the entire 
small F phase). Their effects are expected to round the 
step in Gdiff very near to Fc. 



II. EFFECTIVE MODEL FOR SLIDING CDWS 

Charge-density waves are coupled electron-phonon ex- 
citations which exist in a class of anisotropic metals con- 
sisting of weakly coupled chains. In these materials the 
electronic density is sinusoidally modulated along the 
chain {x~) direction. 



p(r) = po + pi cos(2fcFa; + 4>{r, t)), 



(2.1) 



where kp is the in-chain Fermi wave vector and pi 
the amplitude of the charge modulation. At low tem- 
perature, due to the gap in the dispersion relation for 
the amplitude, amplitude fluctuations are strongly sup- 
pressed and the dynamics can be described in terms of 
the phase (j) only. The Hamiltonian H for a CDW in a 
dimensional metal can be written as |l[] 



n 
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+ drV{r)piC0s{2kFX + <j){r,t)), 



(2.2) 



where we have rescaled coordinates to obtain an isotropic 
elastic term. The coefficient K is a. stiffness constant. 
The effect of impurities is described via a Gaussian ran- 
dom potential V'(r) with zero mean, < ^(r) >= 0, and 
short range correlations, < V{r)V{0) >= Vl3^A(|r|/^o), 



with ^0 a short wavelength cutoff. The overdamped equa- 
tion of motion for the CDW phase variable (p is given by 

= D\/^<j) + u;o + V{r) sm{2kFX + 0(r, t)), (2.3) 

where D = (^)rw?, has the dimension of a diffusion 
constant {l'^t~^) . Here r is the relaxation time of a 
thermally excited phonon, while m and m* are the elec- 
tronic mass and effective mass, respectively. We have 
also let V{r) = {D/ K) V{r ). The second term on the 
right hand side of Eq. ( p. 3D arises from an electric field 
E applied along the a;— direction, with loq = Ee{:^)2kp, 
and can be shifted away hy 4> ^ (p + luqI. The "force" loq 
has dimensions of frequency (t^^) and actually represents 
the "washboard frequency" loq = 2kpv of a freely sliding 
CDW (with velocity v) driven by an external electric field 
E, in the absence of quenched disorder. Equation (2.3) 
is the conventional Fukuyama-Lee-Rice (FLR) model of 
CDW dynamics [0 , which has been studied extensively 
both analytically and numerically, particularly near the 
depinning transition. The FLR model exhibits a depin- 
ning transition at a threshold field Et , corresponding to 
a threshold force ctt = ^T^i:;^)- 

As discussed in Ref. |6j, the FLR equation is incom- 
plete in the strongly driven regime. It is essentially a 
near-equilibrium description, in which only the most rel- 
evant perturbation (the driving field) has been added 
to the equilibrium relaxational dynamics. Several ad- 
ditional effects become important in the sliding state. 

The most important such effect is that of convection. 
In particular, in a CDW mov ing with velocity v, the par- 
tial time derivative dt in Eq. 2.3 must be replaced by the 
total convective derivative Dt = dt+vdx- More generally, 
the linear derivative {dx4>) term arises because the elec- 
tric field breaks the reflection symmetry x — > —x. Note 
that the coefficient of this term is small for small veloci- 
ties, which is why it is neglected in the usual equilibrium 
and near-static (i.e. CDW depinning) contexts. 

A second term ordinarily omitted from the FLR equa- 
tion in equilibrium arises from coupling to the underlying 
periodic lattice. This intrinsic pinning in the direction of 
motion can be incorporated in Eq. ( ^.3[ ) by the replace- 
ment V{r) V{r) + W{x), where W{x) is a periodic 
potential, W{x) = WoCOs{Qx), and Q is in general in- 
commensurate with 2kp. While in other contexts such 
an incommensurate periodic potential can be safely ne- 
glected, we will see that it gives rise to important effects 
for the asymptotic behavior in the strongly driven limit. 

Including both these effects, we arrive at a suitable 
generalization of the FLR equation, 

dt(j) = DV'^ct) + uJo~ adxcj) 

+ [V{r) + W{x)] sm{2kpx + (l){r,t)), (2.4) 

where again a numerical factor has been absorbed into 
the periodic potential W{x). We have allowed for renor- 
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malizations of the convective term by keeping the coef- 
ficient a general, but we expect a '-^ v. Eq. 2A is ca- 
pable of describing the behavior of the CDW (up to the 
aforementioned caveats respecting phase slips and ther- 
mal fluctuations) in the full range of applied fields from 
well below to far above the nominal threshold field. 



In fact, Eq. 2.4 is so general that it is a rather inap- 
propriate point from which to study the moving state. 
This is made evident by making the transformation ^ = 
in order to focus on the fluctuations around the 
uniformly sliding CDW. The resulting equation of mo- 
tion for <\) contains force terms which oscillate rapidly in 
time. To determine their effect at time scales longer than 
2it/ujq, one must develop instead an effective equation of 
motion for a coarse-grained (temporally and spatially av- 
eraged) phase (j). In what follows we will drop the overbar 
and denote the coarse-grained phase variations simply by 

The coarse-graining procedure may be explicitly per- 
formed in two different ways. The simplest method is a 
variant of the high-velocity expansion about the sliding 
state 1^], obtained by iterating a formulation of Eq. 2.4 as 
an integral equation. A more complicated, but conceptu- 
ally more clear approach is to coarse-grain using renor- 
malization group (RG) methods, in which short wave- 
length and high frequency components of cj) are explicitly 
integrated out in a field-theoretic formulation. A similar 
calculation was carried out recently in a different context 
by Rost and Spohn Both approaches are straight- 

forward but tedious, and we simply quote the results in 
what follows H]. 

Several simplifications are obtained in this effec- 
tive coarse-grained description. The most impor- 
tant is the modification of the random potential term 
V{Y)sm{2kpx + (j)), which, as mentioned before be- 
comes oscillatory. A careful treatment reveals, how- 
ever, that this term does not strictly average to zero 
in the coarse-grained model. Instead, as discussed in 
Rcfs. 0,§, it generates an effective spatially varying 
drag force Fp{r). To leading order in — , its correlations 

are < Fp{r)Fp{0) >= FgS{r) with Fq = This may 
be understood physically as simply reflecting variations 
of the impurity density in different regions of the sample, 
which then exert a spatially random drag force on the 
CDW. 

An important difference between this term and the 
original sine-Gordon type term is that it does not pre- 
fer any particular value of the phase variation (p. This 
is in fact an exact result in the moving phase, reflecting 
the non-trivial transformation property (j) + '-^ot un- 
der a time-translation t ^ t + t. In general terms, the 
equilibrium ordered phase of the CDW is described as a 
state of spontaneously broken spatial translation symme- 
try. This state is highly susceptible to disorder, because 
randomness explicitly breaks precisely this symmetry - 
i.e. it acts as a random field. By contrast, the sliding 
CDW breaks time-translation symmetry, which is an ex- 



act invariance of the system, even with V ^Q. 

A second simplification occurs in the intrinsic pinning 
term. Like the random potential, this term also becomes 
oscillatory in time, but generates a non-trivial correc- 
tion upon coarse-graining. To second order in a gradient 
expansion, the correction has the form of an additional 

drag force 5Fw ~ ~ ci|V_l^P — C2|i9a;0P], where 

~ C2 l/(2fci?)^ are constants. Physically, the gradi- 
ent corrections arise because the drag force from intrinsic 
pinning becomes less effective as the CDW wavevector 
(whose local shift is proportional to Vip) becomes less 
commensurate with the underlying lattice. For simplic- 
ity, we will focus on the isotropic case ci = C2, which is 
expected to be approximately correct for CDWs whose 
density profile is well approximated by the single Fourier 
harmonic form of Eq. 



2.1 and which is not too far from 



commensurability. The resulting gradient-squared cor- 
rection is a realization of the Kardar-Parisi-Zhang (KPZ) 
nonlinearity in the CDW system. 

The final coarse-grained equation of motion is 



d,6 = DV-^ 



ad,~^ + Fp{v) + ^{V~^f, (2.5) 



where An 



This coefficient is positive, because 

a mis-oriented CDW (with V0 7^ 0) is less slowed down 
than an aligned one (with V(t> — 0). 

Eq. 2.5 is the basis for our study of the moving state. 



We caution, however, that some information is lost in this 
approach, and various non-universal high energy features 
of the CDW behavior are no longer easily calculable. An 
important example is the full form of the I-V curve. As 
can be explicitly seen in the coarse-graining procedure, 
the CDW frequency dt(j) as a function of E or loq has 
non-trivial contributions from the sh ort- wavelength de- 
grees of freedom not contained in Eq. |2.5| . An additional 
difficulty is that the drag forces Fq and Aq are strongly 
force dependent. Our long- wavelength description does, 
however, capture the singular part of the CDW velocity. 
We define 



SVsing ^ < dt<i 



>, 



(2.6) 



where the overbar denotes a spatial average and the 
brackets denote the disorder average. The quantity 



SVs 



is actually a frequency shift. Note that, because 



this includes only the singular part of the v{ujo) relation, 
there is no particular preferred sign for dSvsmg/duJo- 

The spatial fluctuations of the phase can be charac- 
terized by their growth with the system size. A useful 
measure of such fluctuations employed in the study of 
interface dynamics is the "interface width" in the long- 
time saturated regime, given by 



w{L^,L^) =< [c/.(r, t) - (/.(r, t)]2 >i/2 (2.7) 
in a d-dimcnsional system of size L^L"^^. 
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If the KPZ term is neglected in Eq. (2.5), the equation 
is hnear and can be solved exactly by Fourier transfor- 
mation, as discussed in Ref. The CDW response is 
linear and < dt4> >= ^Oj i-e., Svging—O. The random 
mobility yields a static distortion of the CDW, 



Dq^ + iaqx 



27r<5(w) 



(2.8) 



function 

^ n 



IS 



The corresponding correlation 
|^(q,w)|2 S{q)2TT6{uj), with S{q) = 
static structure function. The (7dx4> term in Eq. ( |2.5| ) 



< 
the 



is crucial in determining the decay of spatial correlation 
in the moving state. If this term is absent, fluctuations 
are isotropic, with S{q) ~ g""*, so that w{L) ~ L^^-d-)/'^^ 
for Lx = L±_ = L. In particular in c? = 1, w{L) ^ L^^^ 
and the system will develop a "groove" instability of the 
type discussed in The case d = 2 is marginal with 
w{L) ~ L. The adxip term suppresses the growth of fluc- 
tuations in the a;— direction. When this term is present, 
in the limit where Lx » L± and cr >> 1, the CDW is 



"riding over" the static disorder and w{Lx,Lji 



-1/2 



III. TILTED MAGNETIC FLUX LINE ANALOGY 



When the KPZ coupling Aq is nonzero, Eq. (2.5) can be 
mapped into the problem of a directed path in a random 
potential via the well known Cole-Hopf transformation 
iQ. By letting ^'(r,t) = e^'^^'^'*\ a linear equation of 
motion for ^'(r, i) is obtained. 



(3.1) 



The solution of Eq. (3.1) can be written as a path- 
integral, 

r TP /dt'[i(<ir/<it')^-o-£-AoFp(r)] 

«'(r,<)= / V[r]e « . (3.2) 

(0,0) 



Eq. (3_^) can also be interpreted as the partition function 
of a tilted magnetic flux line in the presence of colum- 
nar pinning centers. A single magnetic flux line in a 
(d + 1)— dimensional sample of thickness L in the direc- 
tion of the applied field H, chosen as the z direction 
(H — Hqz), is parametrized by its trajectory r(z) as 
it traverses the sample along the field direction. The 
sample contains columnar pinning centers aligned with 
the z— direction that can pin the flux line over its entire 
length. An additional magnetic field Hj^ applied perpen- 
dicular to the z— direction tilts the flux line away from 
the direction of the columnar pins. The fiux line free 
energy is then given by. 



G 



dz 



rdr 

dz 



h 



C/(r(z))}, 



(3.3) 



where ei = ( ^^^ ^ Yln{n) is the tilt modulus (we as- 
sume for simplicity an isotropic superconductor), h = 
H_L(/)o/47r, and U{v) is the random pinning potential gen- 
erated by the columnar defects. The pinning potential is 
correlated along the direction and has short range corre- 
lations in the plane, with < U{y)U{y') >— A6{r — r'). 
The partition function of a vortex line with fixed end 
points r(0) = and r(L) = is obtained by summing 
the Boltzmann factor e-a/ksT ^^^^ ^^^1 paths connecting 
the end points and it is given by 



(r,t) 



z 



V[r]e 



+ U{riz)) 



(3.4) 



(0,0) 



The dynamics of a d— dimensional driven CDW at T = 
can therefore be mapped onto the tilt response of a mag- 
netic fiux line in a (c? -t- 1)— dimensional superconductor 
with columnar pins, at finite temperature. In this map- 
ping the time argument of the CDW corresponds to the 
flux line coordinate z along the field direction, the diffu- 
sion constant D plays the role of temperature, according 
to D ^ ^ , and the driving force a corresponds to a tilt 
field h = hx with a ^ j^. The correspondence between 
the various CDW and flux line quantities is summarized 
in table 1. 



CDW 


Flux line 


D 


r/2e~i 


a 


h/h 


Ao^o 




— Apl^Wsmg 


g/h 


■I dSVsing 

^0 da 





Table 1 



Eqs. (3^) and (3^) differ by a constant term in 
the flux line free energy which represents the field energy 
associated with the tilt field H^. In the absence of a 
tilt field, the flux line is localized on the strongest colum- 
nar defect. At low temperature the localization length, 
deflned as the radius of the tube to which the flux line 
is conflned, is of the order of the range of the pinning 
potential. Thermal fluctuations increase the localization 
length, but are not sufficient to depin the flux line in 
d — 1,2. A sufficiently strong perpendicular field Hj_ 
will, however, depin the flux line. The response of the 
flux line to the fleld is measured by the average induction 
in the direction B± of the transverse field. We define a 
dimensionless induction b i = ^jr . This is also the mean 

novo 

slope of the tilted flux line. The induction can be written 
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aa o I = -t- ^:n'nl i , 
sionless) total magnetization and g{h) is the Gibbs free 
energy per unit length of the tilted flux line. It has been 
argued that a flux-line array pinned by columnar defects 
exhibits a transverse Meissner effect, with b± — for tilt 
field below a critical value he and b± ^ ioi h > h^- The 
tilt response of a single flux line pinned by one columnar 
defect in (1 + 1)— dimensions can be evaluated analyti- 
cally (see Appendix A and jl^). One finds that in the 
limit Lx —> oo there is a transverse Meissner effect for 
h < he- For a (5— function pin with U{r) = —Uo6{r), we 
find he ^ ^ ^^T^ ■ In the pinned configuration for h < he 



Anm^, where m± 



dg 
dh 



is the (dimen- 



the line is 



the flux line free energy per unit length is g 
so that TO^ = —i^ and h±^ = 0. For h > he 

^ ei ^ — 

depinned and 5 = 0. This gives m± — and 6± = 
The transition from a pinned to a depinned configuration 
is associated with a jump discontinuity in the induction 
or tilt slope b± at he and can therefore be classified as a 
first order phase transition. For a general pinning poten- 
tial we estimate he as the field required to depin the flux 
line, hl/2ei ~ VA, or he ~ (2eiVA)i/2. Similar con- 
clusions were reached by Balents and Simon for the 
tilt response of a single flux line in a random distribution 
of columnar pins in (1 -I- 1)— dimension. Also Hatano 
and Nelson very recently related the depinning of 
a flux line from columnar defects by a transverse field 
to the localization transition of a quantum particle in a 
constant imaginary vector potential. By exploiting this 
mapping they showed that the transverse Meissner effect 
persists in both d = 1,2. The question of whether an 
array of many interacting flux lines pinned by columnar 
defects will also exhibits a sharp transition is, however, 
still open. 

The transverse Meissner effect for the vortex line trans- 
lates into a first order dynamical phase transition of the 
sliding CDW. We recall that the CDW driving force a 
corresponds to the tilt field h and Svsing ^ g- There is 
then a transition at a characteristic dc from a state with 
Svsing 7^ for (7 < (Tc to a state where the external drive 
dominates and washes out the effect of disorder in the 
direction of the drive, yielding Sv, 



sxng 



for a < (Tc- 

The driving force Uc where the transition occurs can be 
estimated from the flux-line analogy using table 1 and 
the estimated he « (2eiVA)"^^^ as CTc ~ y^2XoFo. The 
sliding phase with 6vsing 7^ at small driving forces cor- 
responds to the situation where the flux line is pinned on 
the strongest columnar defects and exhibits a transverse 
Meissner effect. As we will show in the next section, this 
is a disorder-dominated regime for the sliding CDW with 
a "rough" spatial profile of the phase (/>(r, t). We will re- 
fer to this phase as a rough sliding phase. For cr > Cc the 
CDW is in a shding phase with Svsing = 0, correspond- 
ing to a flux line depinned by the tilt field and "riding 
over" the columnar pins. As shown below, this phase 
is characterized by anisotropic spatial fluctuations of the 
phase. Fluctuations are suppressed in the direction of the 



driving force and we will refer to this phase as a "flat" 
phase. 



IV. NUMERICAL RESULTS 



We have integrated numerically Eq. (2^) in both 
d = 1, 2 by discretizing the spatial coordinates, with lat- 
tice spacing equal to the range Rp of the pinning poten- 
tial, chosen as our unit of length. We assume an initially 
flat configuration , 0(r, t = 0) = 0, and follow the dy- 
namics until the system relaxes to a steady state. The 
relaxation time scales as L'^L'^~'^. The average CDW 
properties in the saturated sliding state are evaluated 
by performing both a time average and an average over 
many realizations of the disorder. 
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Fig.l 

Figure 1 shows the frequency shift Svsmg of the CDW 



defined in Eq. (2.6) as a function of the driving force 
a for d = 2. A similar behavior is obtained in d = 1. 
As the system size is increased one observes a tran- 
sition between two sliding phases. For the parameter 
values used here the estimated critical force is CTc = V2, 
a value that agrees quite well with our numerical results. 
At large driving forces disorder is washed out by the ex- 
ternal drive and Svsmg approaches zero as 00. For 
small driving forces both impurities and intrinsic pinning 
are important and yield a large Svsing ■ While the drop of 



SVs 



above a critical force does become sharper as 



increases, the approach to the sharp transition expected 
in the limit Lx,L± — > 00 is rather slow. This can be 
understood by examining the dependence on the system 
size Lx of the free energy of a single flux line pinned by a 
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single columnar defect evaluated in Appendix A. As dis- 
cussed in the previous section, this simple model exhibits 
a first-order depinning transition in the limit ^ oo. 
On the other hand, the finite size corrections to the flux 
line free energy are large in the region h > he, as shown in 
Fig. 5. For h < he the flux line is localized on the colum- 
nar pin and does not "see" the rest of the system. As a 
result, in this region the finite size corrections to the free 
energy vanish exponentially with system size. For h > he 
the flux line is delocalized and samples the entire system. 
In this region the free energy is quite sensitive to the fi- 
nite system size, with g{Lx) — g{oo) ~ The scaling of 
SvsingiLx, L±) in our driven CDW problem - which maps 
onto a tilted flux line pinned by a random distribution of 
many columnar defects - is even slower than obtained in 

1/2 

the single pin model , with Svsing ~ Lx approximately. 
This slow approach to the asymptotic limit can, however, 
be understood by the same physical argument. 




F„=1.0 
X„=1.0 

L^=256 Ly=32 □ 
L =32 L =32 i 

X y 

L =4 L =32 a 



rig.2 



at the transition 



To better display the change of Sv^ing 
we show in Fig. 2 the evolution of the curve "'2'^'' versus 
a with system size for d = 2. In the limit ^ oo 
the derivative will exhibit a jump discontinuity at the 
transition. As the driving force cr is proportional to the 
applied voltage (electric fleld) and dvsmg is proportional 
to the current J, this corresponds to a jump discontinuity 
in the differential conductance Gdiff = Our model 
only allows us to predict the singular shift Svsing and 
therefore the singular behavior of Gdiff, but not the full 
nonlinear form of the I(V) curve. The magnitude of the 

jump discontinuity in is, however, proportional to 

the jump discontinuity in Gdifi- The precise relation and 



the possibility of observing this effect will be discussed in 
the next section. From table 1 we see that corre- 
sponds to the magnetization of the tilted flux line. Fig. 
2 displays then the transverse Meissner effect discussed 
earlier. 

By translating the results obtained in Appendix A for 
the ffux line, we find that if the random force Fp{r) of Eq. 



( 3.1 ) is replaced by a single (5— function pin, the frequency 
dv^ng is given by 



2Ao 



(4.1) 



This form is shown as a dashed line in Fig.l and fits very 
well our data for a < ac- This is easily understood be- 
cause in the region a < ac the fiux line is localized onto 
the strongest pin and its free energy is basically unaf- 
fected by the presence of the other defects. 

The two sliding states of the driven CDW differ quali- 
tatively in the behavior of the spatial fiuctuations of the 
coarsed-grained phase (j3{r,t). For a < ac pinning domi- 
nates the dynamics. The sliding state is rough with large 
spatial fluctuations of the phase both in the directions 
parallel and perpendicular to the external drive a. For 
a > ac the term adx4> washes out the effect of pinning in 
the x-direction, damping out phase fluctuations in this 
direction. In this case, the spatial fluctuations of the 
phase are anisotropic and are suppressed in the direction 
of the external drive. This behavior is shown qualita- 
tively in Fig. 3 that displays contour plots of the CDW 
phase for increasing values of a. 

To quantify the behavior of phase fluctuations in the 
two sliding states, we hav e ex amined the interface width 
'w{Lx,L±) defined in Eq. (2/7). In the isotropic disorder- 
dominated phase for a < ac we expect w ^ Lx ^ L±. 
To understand this, we recall that when a = the 



path-integral solution (3.2) of the CDW problem can 
also be interpreted as the transfer matrix solution of the 
Schrodinger equation for a quantum particle in a random 
potential in d spatial dimensions and imaginary time [p^ . 
The width w{Lx,L±) corresponds to the fluctuations in 
the energy of the quantum particle as a function of sys- 
tem size. For d = I the c^uantum particle is always local- 
ized. The states are exponentially localized and one can 
show [|l^ that the energy fluctuations scale as the system 
size, i.e., w{Lx) ~ Lx- A similar behavior is expected for 
d = 2. In the large-cr phase, we postulate an anisotropic 
scaling ansatz for the interface width. 



wiLx,L^)^L^fiLjL^,), 



(4.2) 



where x and z are two new exponents. For >> L± 
the system looks one-dimensional, extended along the 
x— direction. An approximat e eq uation of motion in this 
regime is obtained from Eq. (2.5) with \/ dx- For large 
a both pinning by impurities and intrinsic pinning which 
yields the KPZ nonlinearity are negligible compared to 
the adx4> term and one can obtain an approximate solu- 
tion of the equation, which yields w^Lx, L±) 



rl/2 
^x 
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The scaling function /(s) in Eq. ( |4.2| ) must therefore 

X — 1/2 

obey /(s) s ^ for s << 1. This result is easily 
understood by exploiting the mapping of driven CDW 
dynamics onto the problem of the tilted magnetic flux 
line. For a >> CTc {h >> he) the flux line is delocal- 
ized and "rides over" many columnar defects. The inter- 
face width w{Lx, L±) corresponds the fluctuations of the 
flux line free energy, which in this limit is determined 
by a sum of independent random energies, yielding a 

1/2 

w{Lx, L±) ~ Lx scaling for » L±. In the opposite 
limit, L±^ >> the system looks one-dimensional along 
the y— direction for d = 2. The a-dxtj) term has no effect 
and fluctuations are dominated by disorder, with the re- 
sult w{Lx, Lj_) ^ L±. This follows again from the expo- 
nential localization of states of a one-dimensional quan- 
tum particle in a random potential. The scaling function 
obeys /(s) ~ s for s >> 1 and w{Lx, L±) ~ L~^~^L±^ for 
>> in d = 2. A scaling collapse of our numer- 
ical results for the interface width in d = 2 is shown in 
Fig. 4. The best collapse is achieved with x = 0.5 ± 0.05 
and z = 0.85 ± 0.05. These results agree well with the 
asymptotic values discussed earlier. 



V. CONCLUSION 

We have studied the dynamics of driven CDWs mov- 
ing through a random medium at zero temperature, at 
driving forces well above the threshold force where the 
depinning transition occurs. The CDW model consid- 
ered incorporates new nonequilibrium terms which are 
important in the strongly driven regime and are gener- 



ally not included in the FLR model. We have found 
that in d = 1,2 the driven CDW exhibits a first order 
phase transition at a critical driving force CTc ^ \/2XqFq. 
For a < ac disorder controls the dynamics yielding a 
rough sliding phase with spatial fluctuations of the phase 
that grow linearly with Lx, indicating that the phase- 
only model breaks down. For a > ac, the driving force 
washes out the effect of disorder in the direction of mo- 



tion. The CDW slides uniformly with Svg 



0. This 



moving phase is highly anisotropic as the external drive 
suppresses the spatial fluctuations of the phase in this 
direction. The CDW remains "rough" in the direction 
transverse to the external drive. By using the Cole-Hopf 
transformation, the problem of CDW dynamics at large 
driving force can be mapped onto the problem of the tilt 
response of a magnetic flux-line pinned by columnar de- 
fects. The dynamical transition of the sliding CDW cor- 
responds to the transverse Meissner effect of the tilted 
flux line. 

Our coarsed-grained model of phase dynamics given 
in Eq. ( |2.5| ) only applies in the strongly driven phase, 
well above the T — threshold field, Et- In the weak 
pinning limit, one can relate Et to the Lee-Rice length 
Lq that represents the typical domain size in a pinned 
CDW. The pinning length Lq is the length where the 
elastic strains induced by disorder are of order one and 



is given by Lq = 



( rS-l ) 



Here c is the average 



spacing between the CDW chains and is the range of 
the disorder potential. We expect ~ c- The threshold 
field is then estimated by balancing the total force on a 
domain of size Lq to the elastic force acting on the same 
domain, with the result, Et = ^^pfcf L'^-^. This corre- 

' " ^ 2ec o 

sponds to a threshold force ctt = ETe{-^)- The first 
order transition is predicted to occur at a critical force 



where to leading order in i, Aq ^ 



w, 



8k' 



^ and Fo 



0- 



(5.1) 

4kF<T ' 

]) self- 
(5.2) 



Recalling that a — Eer/m*, we can solve Eq. (^ 
consistently for a critical field Ec, with the result 

m" ' ' e 

The first order transition will only be observable if £'c > > 

Et. We find EJEt ~ ^i^y^lx:) for d = '2- The 
first order phase transition at Ec may therefore be ob- 
servable in a dirty material (Lq « c) with appreciable 
intrinsic pinning [Wq >> Vq). Using the results of our 
calculation of the tilt response of a flux line pinned by a 
single defect (Appendix A), we estimate the magnitude 
of the jump discontinuity in Svsing at ac as |— ~ 
— ^2Fq/Xq. This corresponds to a jump disconti- 
nuity in Gdiff that can be expressed in terms of micro- 

scopic CDW parameters as AGdiff = (||)(t7)^(^) 
for d — 2. The discontinuity in Gdig is very small when 
the condition Wo > > Vq of observibility of transition is 
satisfied. 
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APPENDIX A: SINGLE FLUX LINE AND ONE 
COLUMNAR DEFECT 

In order to gain some understanding of the dependence 
of our results on the size Lx of the system in the direc- 
tion of the external driving force (or tilt field for the 
magnetic flux line), it is instructive to consider the ac- 
tion of a transverse tilt field on a fiux line pinned by a 
single attractive columnar defect. The partition function 
is given by Eq. (3.4), with U{r) the pinning potential 
due to the single impurity, chosen for simplicity as an 
attractive (5-function, 



U{r) = -UoSir). 



(Al) 



Following Ref. |17[] , Eq. ( |3.4| ) can be thought of as a path 
integral in imaginary time and the partition function can 
be rewritten as a quantum mechanical matrix element. 



Z(ri,0;i) =<ri|e-^^/^|0>, 



(A2) 



where |0 > and < r^l are initial and final states local- 
ized at and r_L , respectively, and the "Hamiltonian" Ti, 
is the operator. 



zei 



(A3) 



The operator TC is nonhermitian as TO{h) = 7i(— /i). To 
find its spectrum we need to solve both the right and left 
eigenvalue problem, defined by, 



and 



(A4) 



(A5) 



where u^{x) and u^{x) are the right and left 
eigenf unction, respectively, normalized according to 
dxu^{x)u^{x) = 1, and E„ are the corresponding 
eigenvalues. The path integral (3.4) can then be ex- 
pressed in terms of the eigenvalues and eigenfunctions 
of the fictitious quantum problem, as 



This is also equivalent to writing the path integral in 
terms of the eigenvalues of a corresponding transfer ma- 
trix. In the limit L —^ oo the smallest eigenvalue dom- 
inates and the free energy per unit length g{h) of the 
flux line is determined by the real part of the ground 
state energy Eq of the quantum problem, according to 
g{h) = Eo + h^ /2ei. For localized states, the ground state 
wavefunction uo(r) determines the localization length of 
the flux fine |l|]. 

For simplicity we begin by considering a flux line in 
1-1-1 dimensions, with x the direction of the applied tilt 
fleld. The vortex free energy is given by the ground state 
of the non-hermitian "Schrodinger equation", given by 



2ei ^ dx 



hf - Uq5{x) u{x) = Eu{x), (A7) 



to be solved with periodic boundary conditions, u(0) = 
u{Lx). We are considering here the right eigenvalue prob- 
lem and dropping for simplicity of notatio n th e labels on 
the eigenfunction. The solutions of Eq. ( A_7 ) for x 
are given by u±{x) = e^^^^^^hx/T ^ with k — \J ~2liE jT^ . 



The general solution of Eq. (A?) can then be written as 



(A8) 



where A and B are constants to be determined by the 
periodic boundary condition and the condition that the 
wavefunction has a slope discontinuity determined by the 
(5-function, (|f )^ - (^)^ = 2^y(o). The condition 
for a nontrivial solution to exists yields the eigenvalue 
equation, given by, 

cosh(ft,La;/r) — co^inLx) + — sinh(KL2;). (A9) 

In the limit i^, — > oo, there is one localized ground state 
for h < he, with he = iiUo/T. For h > he all states are 
extended. The ground state energy is given by. 



= h<he 



77>00 

^0 



h > he, (AlO) 

The ground state wavefunction is given by 



u^{x) = 
u^ix) - 



T -ihc-h)x 



1 



h < he 
h > he. (All) 



It is exponentially localized for h < he, with localization 



length ^ 



(ha-h) ■ 



If the system is infinitely long in the 



ficld(z— ) direction (L oo), the fiux line free energy per 
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unit length is simply g{h) = Eq + H'^ /2ei. The free energy 
per unit length is shown in Fig. 5 as a function of h (thick 
line). Clearly the magnetization m± is = —h/ei be- 
low he and cancels the applied field, yielding h±_ = 0, 
as required for transverse Meissner effect. For h > he, 
b± ~ h/ii, which is the value in the absence of disorder. 
The induction has a jump discontinuity at he- 

We now discuss the corrections to the above results 

due to a finite system size L^. This will be useful for un- 
derstanding our numerical results for the driven CDW. 
One can study analytically the finite size corrections in 

the limit L^h » 1. Keeping the leading finite size cor- 
rection, the real part of the ground state energy is given 

by 



Eo{Lx) 



hi 
2h 


l + 2e 




h^ 




hT 


2ii 




h 


h^ 


The 




2ei 







h < he 



he 
h 




Fig.5 



h > he, (A12) 



where we have assumed h » he- In the region h < he, 
where the flux line is pinned on the defect, the finite 
size corrections vanish exponentially as with ^ the 
localization length. For h > he the flux line is depinncd 
and samples the whole sample. In this region the asymp- 
totic {h» he) correction to the ground state energy of 
an infinite system vanishes only as The free energy 
g{h, Lx) is shown in Fig. 5 for a few values of system size 

The CDW frequency shift Svsing is simply propotional 
to the negative of g{h) (see table 1). The similarity in 
the finite size dependence of Figs. 1 and 5 is appar- 
ent. Finally, Fig. 6 shows the dependence of the nor- 
malized probability density in the ground state P{x) = 
u^{x)u^(x) on system size Lx- The tilt induces a mis- 
match of order g-{hc-h)Lx/T ^j^g wavefunction at the 
system boundaries. For h « he this mismatch is expo- 
nentially small and the periodic boundary conditions are 
easily accommodated. For h he the boundary condi- 
tions force a large negative gradient in the wave function, 
as apparent from Fig. 6. Physically this corresponds to 
the fact that for /i ~ /ig and finite Lx long sections of 
the flux line are still pinned on a strong columnar pin, 
yielding large finite size effect in the region above the 
depinning transition. 
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Figure Captions 



Fig. 1 The singular part 6vsing of the CDW ve- 
locity as a function of applied force a for various sys- 
tem sizes. The critical force CTc is estimated to be 
(Tc ~ V2Ao7o = \/2 for the set of parameters used in 
the figure. 

Fig. 2 The derivative of the CDW velocity < 
dt(p > with respect to the apphed forces a. This figure 
can also be interpreted as the transverse magnetization 
m(measured in units of j^) versus transverse magnetic 
field for a flux line in columnar defects. 
Fig. 3 Two-dim ensional CDW phase configurations 
4>{x,y.t) — (j){x,y:t) for various driving forces a at long 
time for Ao = 1.0 and Fq ~ 1.0. The estimated CTc is 
(Tc = V^- The contour plots are (a)cr — 0.1, (b)cr = 1.4 
and (c)(7 = 10.0. The relative value of the CDW phase 
(j) is given by the greyscale intensity, with the bright- 



est spots corresponding to the highest (f> and the darkest 
spots the lowest (j). 

Fig. 4 The figure shows the scaling collapse of the 
i nter face width wi^Lr,., L±) according to the ansatz Eq. 
(4.2) for d ~ 2. The parameter values are indicated and 

cr = 50.0. 

Fig. 5 The free energy per unit length g{h) of a 
tilted flux line pinned by a single defect as a function of 
j^. Both the result in the thermodynamic limit ^ oo 
(thick line) and the finite size results are shown. 

Fig. 6 Probability density P{x) = u^{x)u^{x) for 
various system sizes at (a) h = O.Shc and (b) h = 1.2hc- 
In Fig. (a) the probability density decays exponentially 
to zero over a length a; ~ ^ ^ 5 in the limit oo. 
The finite-size correction are determined by the value of 
P{x) in the central fiat region. They are therefore rather 
small even for not too large a value of L^- In Fig. (b) 
- in the limit oo. This value is shown as 



P{x) = 
an horizontal line for each 
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